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Abstract 

Gravitational waves from inspiraling, compact binaries will be searched for 
in the output of the LIGO/VIRGO interferometric network by the method 
of "matched filtering" — i.e., by correlating the noisy output of each interfer- 
ometer with a set of theoretical waveform templates. These search templates 
will be a discrete subset of a continuous, multiparameter family, each of which 
approximates a possible signal. The search might be performed hierarchically, 
with a first pass through the data using a low threshold and a coarsely-spaced, 
few-parameter template set, followed by a second pass on threshold-exceeding 
data segments, with a higher threshold and a more finely spaced template set 
that might have a larger number of parameters. Alternatively, the search 
might involve a single pass through the data using the larger threshold and 
finer template set. This paper extends and generalizes the Sathyaprakash- 
Dhurandhar (S-D) formalism for choosing the discrete, finely-spaced template 
set used in the final (or sole) pass through the data, based on the analysis 
of a single interferometer. The S-D formalism is rephrased in geometric lan- 
guage by introducing a metric on the continuous template space from which 
the discrete template set is drawn. This template metric is used to compute 
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the loss of signal-to-noise ratio and reduction of event rate which result from 
the coarseness of the template grid. Correspondingly, the template spacing 
and total number J\f of templates are expressed, via the metric, as func- 
tions of the reduction in event rate. The theory is developed for a template 
family of arbitrary dimensionality (whereas the original S-D formalism was 
restricted to a single nontrivial dimension). The theory is then applied to a 
simple post^-Newtonian template family with two nontrivial dimensions. For 
this family, the number of templates J\f in the finely-spaced grid is related 
to the spacing-induced fractional loss C of event rate and to the minimum 
mass Mmin of the least massive star in the binaries for which one searches by 
jV~ 2 X 105(0.1/£)(0.2 M0/Mmin)2'^ for the first LIGO interferometers and 
by A/" ~ 8 X 10^^(0.1//:) (0.2 MQ/M^i^f-'^ for advanced LIGO interferometers. 
This is several orders of magnitude greater than one might have expected 
based on Sathyaprakash's discovery of a near degeneracy in the parameter 
space, the discrepancy being due to this paper's lower choice of Mmin and 
more stringent choice of C The computational power V required to process 
the steady stream of incoming data from a single interferometer through the 
closely-spaced set of templates is given in floating-point operations per second 
by P ~ 2 X 10i°(0.1/£)(0.2 Mo/Mmin)^'^ for the first LIGO interferometers 
and by P ~ 3 X C){<d.2 Me/Mmin)^'^ for advanced LIGO interfer- 

ometers. This will be within the capabilities of LIGO-era computers, but a 
hierarchical search may still be desirable to reduce the required computing 
power. 

PACS numbers: 04.80.Nn, 06.20.Dk, 95.30.Sf, 95.75.Pq 
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I. INTRODUCTION 



Compact binary star systems are likely to be an important source of gravitational waves 
for the broadband laser interferometric detectors now under construction [|I| , as they are the 
best understood of the various types of postulated gravity wave sources in the detectable 
frequency band and their waves should carry a large amount of information. Within our 
own galaxy, there are three known neutron star binaries whose orbits will decay completely 
under the influence of gravitational radiation reaction within less than one Hubble time, and 
it is almost certain that there are many more as yet undiscovered. Current estimates of the 
rate of neutron star/neutron star (NS/NS) binary coalescences based on these (very 
few) known systems project an event rate of three per year within a distance of roughly 
200 Mpc; and estimates based on the evolution of progenitor, main-sequence binaries Q 
suggest a distance of as small as roughly 70 Mpc for three events per year. These distances 
correspond to a signal strength which is within the target sensitivities of the LIGO and 
VIRGO interferometers 0,^. However, to find the signals within the noisy LIGO/ VIRGO 
data will require a careful filtering of the interferometer outputs. Because the predicted signal 
strengths lie so close to the level of the noise, it will be necessary to filter the interferometer 
data streams in order to detect the inspiral events against the background of spurious events 
generated by random noise. 

The gravitational waveform generated by an inspiraling compact binary has been cal- 
culated using a combination of post-Newtonian and post-Minkowskian expansions to 
post^-Newtonian order by the consortium of Blanchet, Damour, Iyer, Will, and Wiseman [§, 
and will be calculated to post'^-Newtonian order long before the LIGO and VIRGO inter- 
ferometers come on-line (c. 2000). Because the functional form of the expected signal is so 
well-known, it is an ideal candidate for matched filtering, a venerable and widely known 



technique which is laid out in detail elsewhere ||T0[ and briefly summarized here: 

The matched filtering strategy is to compute a cross-correlation between the interfer- 
ometer output and a template waveform, weighted inversely by the noise spectrum of the 
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detector. The signal-to-noise ratio is defined as the value of the cross-correlation of the tem- 
plate with a particular stretch of data divided by the rms value of the cross-correlation of the 
template with pure detector noise. If the signal-to-noise ratio exceeds a certain threshold, 
which is set primarily to control the rate of false alarms due to fluctuations of the noise, a 
detection is registered. If the functional form of the template is identical to that of the sig- 
nal, the mean signal-to-noise ratio in the presence of a signal is the highest possible for any 
linear data processing technique, which is why matched filtering is also known as optimal 
filtering. 

In practice, however, the template waveforms will differ somewhat from the signals. True 
gravitational-wave signals from inspiraling binaries will be exact solutions to the Einstein 
equations for two bodies of non-negligible mass, while the templates used to search for these 
signals will be, at best, finite-order approximations to the exact solutions. Also, true signals 
will be characterized by many parameters (e.g. the masses of the two objects, their spins, the 
eccentricity and orientation of the orbit...), some of which might be neglected in construction 
of the search templates. Thus, the true signals will lie somewhat outside the submanifold 
formed by the search templates in the full manifold of all possible detector outputs (see 
Fig. 

Apostolatos ||ll| has defined the "fitting factor" FF to quantitatively describe the close- 



ness of the true signals to the template manifold in terms of the reduction of the signal-to- 
noise ratio due to cross-correlating a signal lying outside the manifold with all the templates 
lying inside the manifold. If the fitting factor of a template family is unity, the signal lies 
in the template manifold. If the fitting factor is less than unity, the signal lies outside the 
manifold, and the fitting factor represents the cross-correlation between the signal and the 
template nearest it in the template manifold. 

Even if the signal were to lie within the template manifold, it would not in general 
correspond to any of the actual templates used to search the data. The parameters describing 
the search templates (masses, spins, etc.) can vary continuously throughout a finite range 
of values. The set of templates characterized by the continuously varying parameters is of 



course infinite, so the interferometer output must be cross-correlated witli a finite subset of 
tfie templates whose parameter values vary in discrete steps from one template to the next. 
This subset (the "discrete template family" ) has measure zero on the manifold of the full set 
of possible templates (the "continuous template family"), so the template which most closely 
matches a signal will generally lie in between members of the discrete template family (again, 
see Fig. |l|). The mismatch between the signal and the nearest of the discrete templates will 
cause some reduction in the signal-to-noise ratio and therefore in the observed event rate, as 
some signals which would lie above the threshold if cross-correlated with a perfectly matched 
filter are driven below the threshold by the mismatch. Thus the spacing between members 
of the discrete template family must be chosen so as to render acceptable the loss in event 
rate, without requiring a prohibitive amount of computing power to numerically perform 
the cross-correlations of the data stream with all of the discrete templates. 

The high computational demands of a laser interferometric detector may in fact make it 
desirable to perform a hierarchical search. In a hierarchical search, each stretch of data is 
first filtered by a set of templates which rather sparsely populates the manifold, and stretches 
which fail to exceed a relatively low signal-to-noise threshold are discarded. The surviving 
stretches of data are filtered by a larger set of templates which more densely populates the 
manifold, and are subjected to a higher threshold. The spareseness of the first-pass template 
set insures that most of the data need only be filtered by a small number of templates, while 
the high threshold of the final pass reduces the false alarm rate to an acceptable level. 

Theoretical foundations for choosing the discrete set of templates from the continuous 



family were laid by Sathyaprakash and Dhurandhar for the case of white noise in Ref. [12 



and for (colored) power-recycling interferometer noise in Ref. ||T3[. Both papers used a sim- 
plified (so-called "Newtonian") version of the waveform which can be characterized by a 
single parameter, the binary's "chirp mass" Ai. Recently, Sathyaprakash |T^ began con- 
sideration of an improved, "post-Newtonian" set of templates characterized by two mass 
parameters. He found that, by a judicious choice of the two parameters, the spacing be- 
tween templates can be made constant in both dimensions of the intrinsic parameter space. 



Sathyaprakash's parameters also make it obvious (by producing a very large spacing in one 
of the dimensions) that a two-parameter set of templates can be constructed which, if it 
does not populate the manifold too densely, need not be much more numerous than the 
one-parameter set of templates used in Refs. |T2|,|13|. 



In this paper I shall recast the S-D formalism in geometric language which, I believe, 
simplifies and clarifies the key ideas. I shall also generalize the S-D formalism to an arbitrary 
spectrum of detector noise and to a set of template shapes characterized by more than one 
parameter. This is necessary because, as Apostolatos [|ll[] has shown, no one-parameter set 
of templates can be used to filter a post-Newtonian signal without causing an unacceptably 
large loss of signal-to- noise ratio. 

In one respect, my analysis will be more specialized than that of the S-D formalism. My 
geometric analysis requires that the templates of the discrete set be spaced very finely in 
order that certain analytical approximations may be made, while the numerical methods 
of Sathyaprakash and Dhurandhar are valid even for a large spacing between templates 
(as would be the case in the early stages of a hierarchical search). The small spacing 
approximation is justified on the grounds that at some point, even in a hierarchical search, 
the data must be filtered by many closely spaced templates in order to detect a reasonable 
fraction (of order unity) of the binary inspirals occurring in the universe within range of the 
LIGO/VIRGO network. 

The rest of this paper is organized as follows: In Sec. ||, I develop my generalized, 
geometric variant of the S-D formalism. I then apply this formalism to the general problem 
of detection of gravitational waves from inspiraling binaries, and develop general formulas for 



choosing a discrete template family from a given continuous template family. In Sec. |TT[ I 
detail an example of the use of my formalism, choosing discrete templates from a continuous 
template family which describes nonspinning, circularized binaries to post^-Newtonian order 
in the evolution of the waveform's phase. I also estimate the computing power required for 
a single-pass (non-hierarchical) search using this discrete template family, and compare to 



the previous work of Sathyaprakash fl^. Finally, in Sec. I summarize my results and 
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suggest future directions for research on the choice of discrete search templates. 



II. THEORY OF MISMATCHED FILTERING 



In this section, a geometric, muhiparameter variant of the S-D formahsm is developed. 
Unless otherwise stated, the following conventions and definitions are assumed: 

Following Cutler and Flanagan |l5l, we define the inner product between two functions 
of time a{t) and b{t) (which may be templates or interferometer output) as 



and Sh{f) is the detector's noise spectrum, defined below. 

The interferometer output o(t) consists of noise n(t) plus a signal As(t), where A is 
a dimensionless, time-independent amplitude and s{t) is normalized such that (s|s) = 1. 
Thus, A describes the strength of a signal and s{t) describes its shape. 

Waveform templates are denoted by u{t;fi,X), where A is the vector of "intrinsic" or 
"dynamical" parameters characterizing the template shape and /j, is the vector of "extrinsic" 
or "kinematical" parameters describing the offsets of the waveform's endpoint. Examples 
of intrinsic parameters A* are the masses and spins of the two objects in a compact binary; 
examples of extrinsic parameters /i* are the time of a compact binary's final coalescence to 
and the phase of the waveform at coalescence $o- 

Templates are assumed to be normalized such that (^(^t, A)|^^(^t, A)) = 1 for all /x and 



Expectation values of various quantities over an infinite ensemble of realizations of the 
noise are denoted by E[ ]. 




(2.1) 



Here a(/) is the Fourier transform of a{t), 




(2.2) 



A. 
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The interferometer's strain spectral noise density Sh{f) is the one-sided spectral density, 
defined by 

E[n(/i)r = - /2)^.(/i) (2.3) 

for positive frequencies and undefined for negative frequencies. The noise is assumed to have 
a Gaussian probability distribution. 

Newton's gravitational constant G and the speed of light c are set equal to one. 

A. Formalism 

In developing our formalism, we begin by defining the signal-to-noise ratio. For any 
single template u{t) of unit norm, the cross-correlation with pure noise {n\u) is a random 



variable with mean zero and variance unity (cf. Sec. II. B. of Ref. |T^, wherein it is shown 
that E[(n|a)(n|6)] = {a\b)). The signal-to-noise ratio of a given stretch of data o{t), after 
filtering by u(t), is defined to be 

^"^f, , = {o\u). (2.4) 
rms {n\u) 

This ratio is the statistic which is compared to a predetermined threshold to decide if a 
signal is present. 

If the template u is the same as the signal s, it optimally filters the signal, and the 
corresponding (mean) optimal signal-to-noise ratio is 

E[p] = E[{n + Au\u)] 

= A. (2.5) 

If the template u used to filter the data is not exactly the same as the signal s, the mean 
signal-to-noise ratio is decreased somewhat from its optimal value: 

E[p] =E[{n + As\u)] 

= A{s\u). (2.6) 
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The inner product {s\u), which is bounded between zero and one, is the fraction of the 
optimal E[p] retained in the mismatched filtering case, and as such is a logical measure of 
the effectiveness of the template u in searching for the signal shape s. 

Now suppose that we search for the signal with a family of templates specified by an 
extrinsic parameter vector fj, and an intrinsic parameter vector A. Let us denote the values 
of the parameters of the actual templates by A(fc)). For example, /i^^^ might be the 

value of the time to of coalescence for the kth template in the family, and /i^;.-) might be the 
phase of the kth template waveform at coalescence. 

The search entails computing, via fast Fourier transforms (FFT's), all the inner prod- 
ucts (o|M(/X(fc), A(fc))) for A; = 1, 2, . . . In these numerical computations, the key distinction 
between the extrinsic parameters /j, and the intrinsic parameters A is this: One explores 
the whole range of values of ^i very quickly, automatically, and efficiently for a fixed value 
of A; but one must do these explorations separately for each of the X(k)- In this sense, 
dealing with the extrinsic parameters is far easier and more automatic than dealing with 
the intrinsic ones. 

As an example (for further detail see Sec. 16.2.2 of Schutz [|TB|), for a given stretch 



of data one explores all values of the time of coalescence (to = /U^) of a compact binary 
simultaneously (for fixed values of the other template parameters) via a single FFT. If we 
write the Fourier transform (for notational simplicity) as a continuous integral rather than 
a discrete sum, we get 

{o\uifi(k), \k))) 

(2.7) 

= lo^df e^2-/*o5*(/)^(/. other fM^k),\k)). 
The discrete FFT yields the discrete analog of the function of to as shown above, an array 
of numbers containing the values of the Fourier transform for all values of to. 

Because, for fixed X^k), the extrinsic parameters /x are dealt with so simply and quickly 
in the search, throughout this paper we shall focus primarily on a template family's intrinsic 
parameters A, which govern the shape of the template. Correspondingly, we shall adopt 
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the following quantity as our measure of the effectiveness with which a particular template 
shape — i.e. a particular vector A(fc) of the intrinsic parameters — matches the incoming signal: 

max A(fc))). (2.8) 

Here the maximization is over all continuously varying values of the extrinsic parameters. 
Then the logical measure of the effectiveness of the entire discrete family of templates in 
searching for the signal shape is 

max [ max (s|M(/i, A(fc)))], (2.9) 

k fj, 

which is simply (p.8| ) maximized over all the discrete template shapes. 

In order to focus on the issue of discretization of the template parameters rather than 
on the inadequacy of the continuous template family, let us assume that the signal shape 
s is identical to some template. The discussion of the preceding paragraphs suggests that 
in discussing the discretization of the template parameters we will want to make use of the 
match between two templates £t(/; /x, A) and £t(/; /x + A/x, A + AA) which we will define as 

M(A,AA)= max(M(^,A)|M(^ + A/x,A + AA)). (2.10) 

This quantity, which is known in the theory of hypothesis testing as the ambiguity function, 
is the fraction of the optimal signal-to-noise ratio obtained by using a template with intrinsic 
parameters A to filter a signal identical in shape to a template with intrinsic parameters 
A + AA. 

Using the match ( |2.10| ) it is possible to quantify our intuitive notion of how "close" two 
template shapes are to each other. Since the match clearly has a maximum value of unity 
at A A = 0, we can expand in a power series about A A = to obtain 

1 / d'^M \ 

M(A, AA) ^ 1 + - — — ^TT^ AA^AA^. (2.11) 

^ ' ' 2 V^AA'SAAJ I Ax^=o ^ ' 

This suggests the definition of a metric 
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1 / d'^M \ 

SO that the mismatch 1 — M between two nearby templates is equal to the square of the 
proper distance between them: 

1- M = gijAX'AX^. (2.13) 

Having defined a metric on the intrinsic parameter space, we can now use it to calcu- 
late the spacing of the discrete template family required to retain a given fraction of the 
ideal event rate. Schematically, we can think of the templates as forming a lattice in the 
A^-dimensional intrinsic parameter space whose unit cell is an A^-dimensional hypercube 
with sides of proper length dl. The worst possible case (lowest E[p]) occurs if the point A 
describing the signal is exactly in the middle of one of the hypercubes. If the templates are 
closely spaced, i.e. dl 1, such a signal has a squared proper distance 

g^JAX'AX^ = N{dl/2f (2.14) 

from the templates at the corners of the hypercube. 

We define the minimal match MM to be the match between the signal and the nearest 
templates in this worst possible case, i.e. the fraction of the optimal signal-to-noise ratio 
retained by a discrete template family when the signal falls exactly "in between" the nearest 
templates. This minimal match is the same quantity that Dhurandhar and Sathyaprakash 
in Ref. |jl3| denote as but since it is the central quantity governing template spacing it 



deserves some recognition in the form of its own name. Our choice of name closely parallels 
the term "fitting factor" FF, which Apostolatos introduced in Ref. |]ll| to measure the 
similarity between actual signals and a continuous template family. 

The minimal match, which is chosen by the experimenter based upon what he or she 
considers to be an acceptable loss of ideal event rate, will determine our choice of spacing 
of the discrete template parameters and therefore the number of discrete templates in the 
family. More specifically, the experimenter will choose some desired value of the minimal 
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match MM] and then will achieve this MM by selecting the templates to reside at the 
corners of hypercubes with edge dl given by 

MM =1- N{dl/2f. (2.15) 

The number of templates in the resulting discrete template family will be the proper volume 
of parameter space divided by the proper volume per template dl^ , i.e. 



^= / I ^ (2-16) 
(2y'(1 - MM) /N) 



B. Inspiraling Binaries Detected by LIGO 

The formalism above applies to the detection of any set of signals which have a functional 
form that depends on a set of parameters which varies continuously over some range. We 
now develop a more explicit formula for the metric, given an analytical approximation to 
the LIGO noise curve and a particular class of inspiraling binary signals. 

We approximate the "initial" and "advanced" benchmark LIGO noise curves by the 
following analytical fit to Fig. 7 of Ref . ^ : 

S,{f) = I + 2[1 + (///o)^]}, / > /. ^2.17) 

where /o is the "knee frequency" or frequency at which the interferometer is most sensitive 
(which is determined by the reflectivities of the mirrors and is set by the experimenters to the 
frequency where photon shot noise begins to dominate the spectrum) and 5*0 is a constant 
whose value is not important for our purposes. This spectrum describes photon shot noise 
in the "standard recycling" configuration of the interferometer (second term) superposed on 
thermal noise in the suspension of the test masses (first term), and it approximates seismic 
noise by setting Sh infinite at frequencies below the "seismic-cutoff frequency" fg. 

Throughout the rest of this paper, the "first LIGO noise curve" will refer to ( |2.17| ) with 
fs = 40 Hz and /o = 200 Hz, and the "advanced LIGO noise curve" will refer to ( p. 17] ) with 
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fs = 10 Hz and /o = 70 Hz. These numbers are chosen to closely fit Fig. 7 of Ref. for 
the first LIGO interferometers and for the advanced LIGO benchmark. In this paper, when 
various quantities (such as the number of discrete templates) are given including a scaling 
with /o, this indicates how the quantity changes while /o is varied but fs/ fo is held fixed. 

At this point it is useful to define the moments of the noise curve ( ^2.171 ), following Poisson 
and Will 0, as 



rfc/fo x""?/^ 

nq) = 5,(/o) dx 



fs/fo Sh{x fo) 

dx ; „ 

J{q)^I{q)/I{7). (2.18b) 



dx — (2.18a) 



The upper limit of integration denotes the coalescence frequency or high-frequency cutoff 
of whatever template we are dealing with, which very roughly corresponds to the last stable 
circular orbit of a test particle in a non-spinning black hole's Schwarzschild geometry. 

For both first and advanced LIGO noise curves, the majority of inspiraling binary search 
templates will occupy regions of parameter space for which fc is many times fo- Because we 
will always be dealing with I{q) for q > 0, and because the noise term in the denominator 
of the integrand in Eq. (|2.18|) rises as for / ^ /o, we can simplify later calculations by 
approximating /c = cxd in the definition of the moments. 

To illustrate the metric formalism, we shall use templates based on a somewhat simpli- 
fied version of the post-Newtonian expansion. Since the inner product (|2.1|) has negligible 
contributions from frequencies at which the integrand oscillates rapidly, it is far more im- 
portant to get the phase of u{f) right than it is to get the amplitude dependence. Therefore, 
we adopt templates based on the "restricted" post-Newtonian approximation in which one 
discards all multipolar components except the quadrupole, but keeps fairly accurate track of 
the quadrupole component's phase (for more details see Sees. II. C. and HI. A. of Ref. [0). 
Applying the stationary phase approximation to that quadrupolar waveform, we obtain 



M(/;A,/x) = r^/6exp z 



Z_$„ + 27r/to + ^(/;A) 
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(2.19) 



up to a multiplicative constant which is set by the condition {u\u) = 1 |T8| . 

The function describing the phase evolution in ( |2.19| ), is currently known to post^- 
Newtonian order for the case of two nonspinning point masses in a circular orbit about each 
other as 



*(/;M,r/) 



128 



r/-i(7rM/)-5/3 



1 + f (1 + ¥^)(^M/)2/3 - WninMf) 



+ 10(3058673 ^ 5429 ^ 617 2w ^ .^4/3 
^^^V 1016064 ^ 1008 ' ^ 144 ' )\'^^^^J I 



(2.20) 



(cf. Eq. (3.6) of Ref. |]T3)- Here the mass parameters have been chosen to be M, the total 
mass of the system, and 77, the ratio of the reduced mass to the total mass. 

The actual amplitude ^ of a waveform is proportional to 1/-R, where R is the distance 
to the source; and this lets us find the relation between minimal match and event rate which 
we will need in order to wisely choose the minimal match. Assuming that compact binaries 
are uniformly distributed throughout space on large distance scales, this means that the 
rate of events with a given set of intrinsic parameters and with an amplitude greater than 
A is proportional to IjA'^. Setting a signal-to-noise threshold po is equivalent to setting a 
maximum distance Rq oc 1/^ to which sources with a given set of intrinsic parameters can 
be detected. Thus, if we could search for signals with the entire continuous template family, 
we would expect the observed event rate to scale as I/pq. This ideal event rate is an upper 
limit on what we can expect with a real, discrete template family. 

We can obtain a lower bound on the observed event rate by considering what happens if 
all signals conspire to have parameters lying exactly in between the nearest search templates. 
In this case, all events will have reduced signal-to-noise ratios of MM times the optimal 
signal-to-noise ratio A. This is naively equivalent to optimally filtering with a threshold of 
Po/MM, so a pessimistic guess is 



event rate oc 



V Po 



(2.21) 
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for a fixed rate of false alarms. 

In real life, po is affected by the total number of discrete templates and by the minimal 
match of the discrete template family. This can be seen by the fact that the signal-to-noise 
ratio, 



p= max [ max (o|m(/x, A(A;)))], 

k fj, 



(2.22) 



is the maximum of a number of random variables. The covariance matrix of these variables 
will be determined by the minimal match, and will itself determine the probability distri- 
bution of p (in the absence of a signal) which is used to set the threshold po in order to 
keep the false alarm rate below a certain level. However, since these effects are fairly small 
at high signal-to- noise ratios (such as those to be used by LIGO) and the issue of choosing 
thresholds is a problem worthy of its own paper |jl9[ , we will use (|2.21| ) for the rest of this 
paper. 

For this two-parameter template family, the formula for the match ( p.lO| ) can be simplified 
somewhat by explicitly performing the maximization over the extrinsic parameters /x and 



Afi. Since the integrand in the inner product (2J.) depends on ^t and Afx as exp i[2-KfAto — 
A$o]; there is no dependence on /j, but only on Afi. Maximizing over A$o is easy: instead 
of taking the real part of the integral in the inner product ( |2.1| ), we take the absolute value. 

To maximize over Ato we go back a step. Let us define A*^ = to, and consider the 
(A^ + l)-dimensional space formed by A° and A-'. We expand the inner product between 
adjacent templates to quadratic order in A A" to get a preliminary metric 'jap, where the 
Greek indices range from to (Latin indices range from 1 to A^): 



Q2 



lo^df^-^exp ^[27r/Ato + AvI/(/;A^AA^)] 



-7/3 



OAX^dAX/^ 



-7/3 



ShU) 



f 



AA" = 0' 
(2.23) 



Here A$ = $(/; A^' + AA-') - ^(/; A-') 

We define the moment functional such that, for a function a, 
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J\a\ 



rfc/fo 



. , fs/h Sh{x fo) 



X 



-7/3 



-a{x), 



(2.24) 



and thus 



J 



E 



(Xyi Ob 



^anJ{^ - 3n). 



(2.25) 



We also define the quantities ifja such that 

V'o = 2vr/, 



(2.26) 



where the derivative is evaluated at AA-' = and \E' is the part of ^ in Eq. ( ^.19| ) that is 
frequency dependent (any non-frequency-dependent, additive parts of \E' are removed when 
we take the absolute value in the maximized inner product). Evaluation of the derivative in 
Eq. ( p. 23 ) then shows that, in the limit /c//o oo. 



(2.27) 



Finally, we minimize 7a/3AA"AA^ with respect to Ato (i-e., we project 7a/3 onto the 
subspace orthogonal to to) and thereby obtain the following expression for the metric of our 
continuous template family: 



Qij — lij 



ToiToj 
7oo 



(2.28) 



By taking the square root of the determinant of this metric and plugging it into Eq. ( [^.16[ ), 
we can compute the number of templates M that we need in our discrete family as a function 
of our desired minimal match MM, or equivalently of the loss of ideal event rate. 



III. EXAMPLE: CIRCULARIZED NONSPINNING BINARIES TO 
POST^-NEWTONIAN ORDER 

Although the phase of the inspiraling binary signal has recently been calculated to post^- 
Newtonian order ||^, it is useful to calculate the number of templates that would be required 
in a universe where the waveforms evolve only to post ^-Newtonian order and all binaries 
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are composed of nonspinning objects in circular orbits. There are several reasons for this 
exercise. 

1. Apostolatos [|Tl|] has shown that amplitude modulation of the waveform due to spin 
effects is important in an inspiraling binary search only for a few extremal combinations 
of parameters, and also that (at higher post-Newtonian order) templates without spin- 
related phase modulation can match phase modulated signals almost as well as can 
templates that include spin parameters. Therefore the bulk of the final set of templates 
actually used when the detectors come on-line will not need to include the extra spin 
parameters, and we may ignore them in this preliminary work. 

2. We assume circular orbits because gravitational radiation reaction circularizes most 
eccentric orbits on a timescale much less than the lifetime of a compact binary ||2T| . 



3. The phase of the templates is truncated at post^-Newtonian order for simplicity. Al- 
though Apostolatos has demonstrated in Ref. |Tl| that post^-Newtonian templates will 



not have a large enough fitting factor to be useful, consideration of such a set is a first 
step toward obtaining an adequate set of templates — and it is a particularly impor- 
tant step since the metric coefficients will turn out to be constant over the template 
manifold. 

A. Calculation of the 2-Dimensional Metric 

Having chosen as the continuous template family the set of post ^-Newtonian, circular, 
spinless binary waveforms, we must now choose the discrete templates from within this 
continuous family. The first step is to calculate the coefficients of the metric on the two- 
dimensional dynamical parameter space. 

It is convenient to change the mass parameterization from the variables (M, 77) to the 
Sathyaprakash variables [|1^ 
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ri = :^r'M'''^nh)-"\ (3.1) 

Note that ri and T2 are simply the Newtonian and post^-Newtonian contributions to the 
time it takes for the carrier gravitational wave frequency to evolve from /o to infinity. The 
advantage of these variables is that the metric coefficients in (ri, T2) coordinates are constant 
(in the limit fc ^ /o) for all templates. This is because the phase of the waveform u{f) is 
linear in the Sathyaprakash variables, and so the integral in the definition of the match ( p.lO| ) 
depends only on the displacement (Ari,Ar2) between the templates, not on the location 
(ti,t2) of the templates in the dynamical parameter space. 

The dynamical parameter-dependent part of the templates' phase is given by [Eq. ( ^.20| ) 



truncated to first post-Newtonian order and reexpressed in terms of (ri, T2) using Eqs. ( pJ] 



and (|33)] 



^(/; ri, rs) = ^7r/o(///o)-^/Vi + 2nfoif/fo)-'r2, (3.3) 



and it is easy to read off ipi and 1IJ2 [Eq. (|2.26|) j as the coefficients of ri and T2 |2^. By 



inserting these ipj into Eq. ( p.25|) , the relevant moment functionals can be expressed in 



terms of the moments of the noise: 





= 27r/o J(4), 


jm 


= 27r/o |J(12), 




= 27r/o J(10), 




= (27r/o)2 J(l), 




= (2vr/o)2 f J(9), 




= (27r/o)2 J(7), 




= (27r/o)2 |J(17), 




= (27r/o)^ |J(15), 




= (27r/o)2 J(13). 



(3.4) 



We can compute the needed moments of the noise by numerically evaluating the integrals 
(|]T|). By setting the upper limit of integration to infinity, i.e. by approximating /c//o as 



infinite for all templates under consideration, we find that the moments have the constant 
values given in Table |I|; and therefore the moment functionals ( |3.4| ) have the constant values 
given in Table |I|. Inserting these values into Eqs. ( |2.27| ) and ( |2.28| ) yields, for the coordinates 
(A" = to, = Ti, = T2), the 3-metric and 2- metric 



^ -F0.208 -0.220 -0.168^ 



7«;3 = (27r/o) 



+0.784 +0.481 
+0.309 



and 



(3.5) 



gij = (2vr/o 



0.552 0.304^ 



0.173 



for the first LIGO noise curve, and 



+0.209 -0.257 -0.183^ 



'Jal3 = (27r/o) 



+ 1.320 +0.712 
+0.407 



and 



(3.6) 



(3.7) 



gij = (2vr/o) 



^ 1.01 0.486^ 



V 



0.246 



(3.8) 



for the advanced LIGO noise curve (where the dots denote terms obtained by symmetry). 

We shall also estimate the errors in the metric coefficients due to the approximation 
fc/fo ^ 00: The moment integrals defined in Eq. ( p.l8| ) can be rewritten as 

5x~^/^ .00 5x~^/^ 



dx- 



dx- 



//,//o x-^ + 2(1 + x^) Jf^/fo x-^ + 2(1 + x2) ' 

where the first integral is the expression used in the above metric coefficients and the second 
is the correction due to finite fc/fo- The second integral can be expanded to lowest order in 
fo/fc as 

5 



2(l + g/3) 



ifo/fc 



a+g/3 
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and from this the errors in the moments (and therefore in the metric coefficients) due to 
approximating fc as infinite are estimated to be less than or of order ten percent for the 
first LIGO interferometers and one percent for the advanced LIGO interferometers over 
most of the relevant volume of parameter space. Since the two-parameter, post ^-Newtonian 
continuous template family is known to be inadequate for the task of searching for real 
binaries, these errors are small enough to justify our use of the fc^oo approximation in 
this exploratory analysis. 



B. Number of Search Templates Required 

Since the metric coefficients are constant in this analysis, the formula for the required 
number of templates [Eq. (|2.16| )] reduces to 



U ^ 4^^^ I dr, dr,. (3.9) 



2(1 -MM) 

The square root of the determinant of the metric is given by (27r/o)^ 0.108 for the advanced 
LIGO noise curve and by (27r/o)^ 0.058 for the initial LIGO noise curve, so once we have 
decided on the range of parameters we deem astrophysically reasonable we will have a 
formula for A/" as a function of MM. 

The most straightforward belief to cherish about neutron stars is that they all come 
with masses greater than a certain minimum Mmin, which might be set to 0.2 Mq (based on 
the minimum mass that any neutron star can have [Q) or 1.0 Mq (based on the observed 



masses of neutron stars in binary pulsar systems [2^|). In terms of the variables (M, 77) the 
constraint Mi > M^in and M2 > M^i^ is easily expressed as 

1 



2-M(l - ^1 - At]) > M^in, 

but in terms of the Sathyaprakash variables [Eqs. ( ^.Ij ) and ( |3.2| )] the expression becomes 
rather unwieldy to write down. However, see Fig. ^for a plot of the allowed region in (ri, T2) 
coordinates. 
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For this reason we have found it convenient to use a Monte Carlo integration routine 125 



to evaluate the coordinate volume integral / dri dT2. The Monte Carlo approach becomes 



especially attractive when evaluating the proper volume integral / dri dT2 ydet \\gij\\ for 
cases where the integrand is allowed to vary — and in fact may itself have to be evaluated 
numerically, as will be the case for a post^-Newtonian set of templates. The integral has 
numerical values of 0.18 and 24 seconds^ for initial and advanced LIGO interferometer 
parameters, respectively, assuming a Mmin of 0.2 Mq and arbitrarily large Mmax- The 
integral can be shown (numerically) to scale roughly as /q"^'^ (independent of /o) and as 
^min for ^min ranging from 0.2 to 1.0 solar masses (the dependence on Mmax is negligible 
for any value greater than a few solar masses). 

Inserting the above numbers into Eq. ( p.9|) , we find that 

V 0.03 / V 0.2 M© / ^ ' 



for the first LIGO noise curve and 



AT ^ 8.4 X lOM [JI^I^] (3.11) 

V0.03/ 

for the advanced LIGO noise curve. The fiducial value of MM has been chosen as 0.97 to 
correspond to an event rate of roughly 90 percent of the ideal event rate [cf. Eq. (|2.21| )]. 
In terms of the template-spacing-induced fractional loss C of event rate, the number of 
templates required is 

AT ~ 2.4 X 10^ f — V7^^) " (3.12) 



0.1/ \0.2Mfr 



for the first LIGO noise curve and 



0. 



for the advanced LIGO noise curve. 
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C. Template Spacing 



With the aid of the metric coefficients given in Eqs. ( |3.(i| ) and ( ^.81 ), it is a simple task 
to select the locations of the templates and the spacing between them. 

Because the metric coefficients form a constant 2x2 matrix, we can easily find the 
eigenvectors e^.^ and e^.^ of \\gij\\ and use them clS cLXGS to lay out a grid of templates. The 
numerical values are 

e^, = 0.874e^, + 0.485e„, 

(3.14) 

= -0.485e^j + 0.874e^2 
for the ffist LIGO noise curve and 

e^, = 0.899e^, + 0.437e„, 

(3.15) 

= -0.437e^, + 0.8996^^ 
for the advanced LIGO noise curve. The infinitesimal proper distance is given in terms of 
the eigen-coordinates as Ei{dxiy + E2{dx2)'^, where Ei and E2 are the eigenvalues of the 
metric. 

Therefore we simply use Eq. (p.l5| ) to obtain the template spacings 



2(1 -MM) 

dx, = . ^— ^, J = 1, 2. (3.16) 

We find that the eigenvalues of the metrics (|3.6| ) and (|3.8|) are (27r/o)^ times 0.721 and 
0.00427 (ffist LIGO), and (27r/o)2 times 1.25 and 0.00984 (advanced LIGO). Therefore the 
template spacings are given by 



dxi = 0.22 ms 



0.03 / V 200 Hz 



-1 



,,, = 2.9ms(i^)"7^) (3.17) 
^ 0.03 y 1 200 Hz/ ^ ^ 



for the ffist LIGO noise curve and by 



0.49 msfi^) "7-4-' 

V 0.03 / V 70 Hz 



rfa;2 = 5.6msf^ — MM\ / ( (3.18) 
V 0.03 / V^OHzy ^ ^ 
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for the advanced LIGO noise curve. Figure |^ shows the locations of some possible templates 
superposed on a contour plot of the match with the template in the center of the graph. 



D. Computing Power Requirements 



Drawing on the previous work of Schutz |16| concerning the mechanics of fast-Fourier 



transforming the data, we can estimate the CPU power required to process the interferometer 
output on-line through a single-pass (non-hierarchical) search involving M templates. 

Although the data will be sampled at a rather high rate (tens of kHz), frequencies 
above some upper limit fu — 4/o can be thrown away (in Fourier transforming the data) 
with only negligible effects on the signal-to-noise ratio. This lowers the effective frequency of 
sampling to 2/„ (the factor of two is needed so that the Nyquist frequency is /„) , and thereby 
considerably reduces the task of performing the inner-product integrals. If the length of the 
array of numbers required to store a template is F and that required to store a given stretch 
of data is D, the number of floating point operations required to process that data stretch 
through M filters is 

Ni,o. ~ DJ\f{lQ + SlogaF) (3.19) 



[cf. Eq. (16.37) of Schutz [|T6l, with the fractional overlap between data segments x chosen 
as roughly 1/15]. 

Actually, F varies from filter to filter, but most of the search templates occupy regions 
of parameter space where the mass is very low — and thus the storage size of the template, 

i^^2/„n(/o//,)«/^ (3.20) 

is very large The longest filter is the one computed for two stars of mass Mmin, so by 
inserting rj = 1/4 and M = 2M^i^ into Eq. (U) and combining with Eq. (|3l20| ), we find 
that we can make a somewhat pessimistic estimate of the computational cost by using 

14 / f f (3 21) 

\0.2 Mq] \^10Hz; \70RzJ ^ ^ 
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The required CPU power V for an on-line search is obtained by dividing Eq. ( |3.19| ) by 
the total duration of the data set, 

to find that 

P~2Ar/„(16 + 31og2F). (3.22) 



Combining Eq. ( |3.22| ) with Eqs. ( ^.101 ) and ( |3.11| ) gives us 



^ \ ^ ( ^min \ / /o ^ 



P ~ 20 Gflops — ^^^^ — ^ (3.23) 
for an on-line search by the first LIGO interferometers and 

^ \ ^ ( ^ra\n \ ( fo ^ 



p. 270 Gflops ^^j (3,24) 

for the advanced LIGO interferometers. 

Although the estimates in the paragraph above are not to be believed beyond a factor 
of order unity, the magnitude of the numbers shows that a hierarchical search strategy 
may be desirable to keep the computing power requirements at a reasonable level for non- 
supercomputing facilities. That is, the data would first be filtered through a more widely 
spaced (low minimal match) set of templates with a relatively low signal-to-noise threshold, 
and only the segments which exceed this preliminary threshold would be analyzed with the 
finely spaced (high minimal match) templates. 

The metric-based formalism of this paper only holds for the finely spaced set of templates 
used in the final stage of the hierarchical search; the template spacing used in the earlier 
stages of the search will have to be chosen using more complex methods such as those of 



Sathyaprakash and Dhurandhar [|T2|,|T3|. 



E. Comparison With Previous Results 



The only previous analysis of the problem of choosing the discrete search templates 
from the two-parameter, restricted post ^-Newtonian continuous template family is that of 
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Sathyaprakash |T^, in which he found that the entire volume of parameter space correspond- 



ing to Mjnin = 1 Mq could bc covered by a set of templates which vary only in T1+T2 — thereby 
reducing the effective dimensionality of the mass parameter space to one. This implied a 
value of M similar to that obtained in the one-parameter (Newtonian template) analysis of 
Dhurandhar and Sathyaprakash in Ref. |T3 . 



It is not possible to fairly compare my value for M to the values given by Dhurandhar and 
Sathyaprakash in Table II of Ref. due to our differing assumptions concerning the sources 



and the desirable level of the minimal match. Therefore I will compare the assumptions. 

Dhurandhar and Sathyaprakash typically consider a minimal match of 0.8 or 0.9 rather 
than 0.97, This would lead to a loss of thirty to fifty percent of the ideal event rate [cf. 
Eq. (|2.21|) ]. If the current "best estimates" of inspiraling binary event rates are correct, 



the ideal event rate for LIGO and VIRGO will not be more than about one hundred per 
year even when operating at the "advanced interferometer" noise levels, and the loss of up 
to half of these events would be unacceptable. 

From Eqs. ( p. 101 ) and ( |3.11D it can be seen that the dependence of N on Mmin is the 



most important factor influencing the number of templates. The two-parameter analysis 
of Sathyaprakash |]T1| uses a value for Mmin of 1 M©, which is based on the statistics of 



(electromagnetically-) known binary pulsars. However, because there is no known, firm 
physical mechanism that prevents neutron stars from forming with masses between 0.2 and 
1 Mq, LIGO and VIRGO should use a discrete template family with Mmin = 0.2 Mq. After 
all, laser interferometer gravitational wave detectors are expected to bring us information 
about astronomical objects as yet unknown. 

During the final stages of completion of this manuscript, a new preprint by Balasubra- 
manian, Sathyaprakash, and Dhurandhar appeared in the xxx.lanl.gov archive |^, applying 
differential geometry to the problem of detecting compact binary inspiral events and extract- 
ing source parameters from them. The preprint applies the tools of differential geometry 
primarily to the problem of parameter measurement rather than that of signal detection, 
and so does not develop the geometrical formalism as far as is done in Sec. |T| of this pa- 
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per. The metric constructed in Ref. |^ is identical to the information matrix which was 
suggested for use in the construction of a closely-spaced discrete template family in the au- 
thors' previous work |T^. While this is quite useful for parameter measurement, it neglects 
maximization over kinematical parameters and thus is not very useful for the construction 
of search templates. Also, the assumptions about MM and Mmin are no different from those 



made in previous work up to and including Ref. Jill, ^"^^ result for M is no different. 

The main difference between the results of Ref. and previous analyses by the same 
authors — and therefore the most important part of the preprint as far as the detection 
problem is concerned — is the introduction of the possibility of choosing search templates 
to lie outside the manifold of the continuous template family. Using an ad hoc example, 
the authors show that such a placement of templates can result in a spacing roughly double 
that between discrete templates chosen from the manifold formed by the continuous template 
family. My analysis in this paper does not consider this possibility, but the formalism of 
Sec. can easily be extended to investigate this problem in the future. 



IV. CONCLUSIONS 

A. Summary of Results 

This paper has presented a method for semi-analytically calculating the number of tem- 
plates required to detect gravity waves from inspiraling binaries with LIGO as a function 
of the fraction of event rate lost due to the discrete spacing of the templates in the binary 
parameter space. This method, based on differential geometry, emphasizes that ultimately 
a finer template spacing is required than has previously been taken as typical in the litera- 
ture, in order to retain a reasonable fraction of the event rate. This paper details the first 
calculation of this kind that uses post-Newtonian templates and a noise curve which takes 
into account the coloration of noise in the detector due to both standard recycling photon 
shot noise and thermal noise in the suspension of the test masses. 
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The result is that it is possible to search the data for binaries containing objects more 
massive than 0.2 Mq thoroughly enough to lose only ~ 10 percent of the ideal event rate 
without requiring a quantum leap in computing technology. The computational cost of such 
a search, conducted on-line using a single pass through the data, is roughly 20 Gflops for the 
first LIGO interferometers (ca. 2000) and 270 Gflops for the advanced LIGO interferometers 
(some years later). This is feasible (or very nearly feasible) even for a present-day super- 
computing facility, but a hierarchical search strategy (using as its first stage a widely-spaced 
set of templates similar to that analyzed by Sathyaprakash [Q) may be desirable to reduce 
the cost. 



B. Future Directions 

A thorough investigation of hierarchical search strategies is in order: How should the 
threshold and the minimal match of the first stage be set in order to minimize the CPU 
power required while keeping the false alarm and false dismissal rates at acceptable levels? 
How would non-Gaussian noise statistics affect the first stage threshold and minimal match? 
Would a hierarchical search benefit by using more than two stages? How is the threshold 
affected by the minimal match when the approximation of high signal-to-noise ratio can no 
longer be made? 

The formalism of this paper should be applied to choose discrete templates from a better 
continuous template family than the one considered here. The best two-parameter templates 
will be based on the highest post-Newtonian order computations that have been performed 
for circularized, spinless binaries, augmented perhaps by terms of still higher order from the 
theory of perturbations of Schwarzschild or Kerr spacetime. I plan to soon apply my geo- 
metric formalism to the post^'^-Newtonian templates which are currently the best available. 
The areas of parameter space where spins cannot be neglected (noted by Apostolatos in 



Ref. |TT|) must also be investigated, and the inclusion of an orbital eccentricity parameter 



should be considered. 
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There are several more issues which I plan to address using my formalism or some 
extension of it. An analysis needs to be made for the case when the signal is not identical to 
some member of the continuous template family (i.e. the fitting factor is not equal to one); 
and the result of such an analysis should be used to set definite goals for both the fitting factor 
and the minimal match in terms of event rate. The effect of non-quadrupolar harmonics of 
the gravitational wave on the construction of search templates should be considered. These 
harmonics have been ignored in all previous analyses of detection and even of parameter 
measurement, but they may have a noticeable effect when a very high minimal match is 
desired. Finally, a systematic investigation of the optimal choice of search templates outside 
the continuous template family is in order. This problem has been addressed in a preliminary 
way in Ref. p^, but is deserving of further scrutiny. 
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FIGURES 

FIG. 1. A schematic depiction of the manifold formed by the continuous template family, here 

represented as a two-dimensional surface lying within a three-dimensional space. The discrete 
template family, shown by the dots, resides within this manifold. The x indicates the location of 
an actual signal, which because it is an exact solution to the Einstein equations does not lie within 
the manifold. The -|- marks the spot in the manifold which is closest to (has the highest inner 
product with) the signal. In general, this location falls in between the actual discrete templates. 

FIG. 2. The two-dimensional region of parameter space inhabited by binaries composed of 
objects more massive than 1 Mq. ti and T2 are expressed (in milliseconds) for /o = 200 Hz, but 
the shape of the region does not change with /q. The upper boundary of the wedge is set by Mmin- 
The left-hand boundary is set by Mmax, but is essentially identical to the T2 axis for M^ax greater 
than a few solar masses. The region below the wedge corresponds to > 1/4, which is a priori 
impossible. 

FIG. 3. Locations of various discrete templates for the first LIGO interferometers are shown by 
dots. The contours indicate the value of the match between a template located at (ri,T2) and the 
template located in the center of the figure. The contours are drawn at match values of 0.97, 0.98, 
and 0.99. The semimajor and semiminor axes of the contour ellipse (along which the templates 
are placed) do not appear to be perpendicular because of the aspect ratio used to make the graph 
readable. 
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TABLES 

TABLE I. Numerical values of the moments of the noise in the limit /c//o — ^ oo, for the 

noise curves of the first LIGO interferometers (fs/fo = 1/5) and advanced LIGO interferometers 



{fs/fo = 1/7). 


Noise moment 


Tr^-vci+ Tir^r^ Tmin/-* 
r irst ijiicrU value 


Advanced LIGO value 


J(l) 


1.27 


1 o£; 
i.ZD 


J (4) 


(J. 927 


U.919 




1 (exact) 


1 (exact) 


J(9) 


1.24 


1.26 


J(10) 


1.44 


1.49 


J(12) 


2.13 


2.31 


J(13) 


2.69 


3.03 


J(15) 


4.67 


5.80 


J(17) 


8.88 


12.7 
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TABLE II. Numerical values of the moment functionals, under the same assumptions as in 

Tal)le I. 



Moment functional 


First LIGO value 


Advanced LIGO value 


J[M 


(27r/o)0.927 


(27r/o)0.919 




(27r/o)1.28 


(27r/o)1.38 




(27r/o)1.44 


(27r/o)1.49 




(27r/o)2l.27 


(27r/o)2l.26 




(27r/o)20.743 


(27r/o)20.756 




(27r/o)2 (exact) 


(27r/o)2 (exact) 




(27r/o)23.20 


(27r/o)24.56 




(27r/o)22.80 


(27r/o)23.48 


jm 


(27r/o)22.69 


(27r/o)23.03 
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